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Abstract 

Using a three stage model of hadron formation we calculate the change of the trans- 
verse momentum distribution of hadrons produced in semi-inclusive deep inelastic 
scattering (SIDIS) on nuclei. In the first stage after its interaction with the virtual 
photon, the struck quark propagates quasi free in the nuclear environment under- 
going multiple collisions with nucleons. During this stage it can acquire transverse 
momentum. In the second stage a prehadron is formed which has a very small elas- 
tic cross section with the nucleons. In the third stage the prehadron turns into a 
hadron. For HERMES energies, prehadron elastic scatterings contribute little to p_L- 
broadening. The acquired extra Ap^ of hadrons can therefore be deduced entirely 
from the first stage of quasi free quark propagation with quark- nucleon collisions. 
We use this model to describe 7r-production on Ne, Kr, Xe and compare with the 
most recent HERMES preliminary data. 



1 Introduction 



Fragmentation of quarks and gluons into hadrons is a consequence of color con- 
finement and is therefore one of the most-interesting parts of non-perturbative 
QCD. Use of nuclear targets allows the experimentalist to position detectors 
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(nucleoli targets) of the hadronization process near the interaction point and 
probe hadron formation on length scales of a few Fermi. 

A struck quark originated from deep-inelastic scattering (DIS) propagates 
through the nucleus experiencing multiple interactions which cause induced 
energy loss. In the string model the medium- induced energy loss results from 
production of new strings in multiple collisions of the quark [1]. In a perturba- 
tive description the excess of energy loss is related to medium- induced gluon 
radiation resulting from broadening of the quark transverse momentum [2] . In 
both approaches induced energy loss rises quadratically with the path-length 
of the quark. 

Although induced energy loss contributes to the attenuation of the produced 
hadron, even stronger suppression may result from absorption of the produced 
hadron, if the struck quark is neutralized and a colorless dipole is formed inside 
the nucleus. In the string model, the production of a hadron may be described 
as a two stage process [3]l^ll5] . First a prehadron (e.g. a dipole of small size) 
is created, which has a reduced absorption cross section, and later the hadron 
wave function is formed. These models [3|^f5] put the main emphasis on the 
effects of absorption and neglect the induced energy loss. The attenuation 
of the produced hadron can be described more consistently within the color- 
dipole approach employing a path integral technique [6|7] which describes the 
evolution of the dipole propagating through the medium. The induced energy 
loss evaluated perturbatively in Ref. [2] was included in Ref. [7]. 

Some models put the main emphasis on the effects of induced energy loss, 
assuming that color neutralization of the quark always occurs at long dis- 
tance from the DIS location, i.e. outside of the nucleus pi9] . Of course, this 
assumption may only be valid in certain kinematical domains. 

A novel mechanism of hadron attenuation was found recently in Ref. [10] . It 
turns out that even if no induced energy loss occurs, and color neutralization 
happens outside of the nucleus, a significant nuclear suppression can result 
from quantum-mechanical interferences between different amplitudes. 

In this paper we study the recent HERMES preliminary results [11] which give 
preliminary data on p^-broadening of pions produced in deep inelastic lepton- 
nucleus scattering (DIS) on Ne, Kr and Xe nuclei. Transverse broadening 
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Ap^ = {pj^A — {p'i)d in the nuclear medium of a nucleus with mass number 
A has been measured as a function of the hadron fractional momentum Zh, 
of the virtual photon energy u and its virtuality Q^, for different nuclear size 
and for different hadrons. Here, the index D in the definition of the transverse 
broadening refers to a deuteron target. 

We use the absorption model [3|5] to calculate the nuclear modifications of 
hadron production in DIS. The hadron formation time is computed analyti- 
cally in the framework of the LUND string fragmentation model as a three-step 
process. In the first stage the quark (or antiquark) ejected from the nucleon 
propagates and undergoes multiple collisions in the nucleus. In the second 
stage color neutralization takes place and a prehadron is formed. Inelastic in- 
teractions of the prehadron or hadron result in a considerable shift of the final 
(detected) hadron towards smaller Zh- We treat this process as absorption. 
Only elastic rescatterings preserve the (pre-)hadron and contribute to broad- 
ening. However, the elastic cross sections of hadrons and prehadrons are very 
small compared to the inelastic cross sections. In the third stage the final state 
hadron is formed from the surviving prehadrons. With an inelastic prehadron 
nucleon cross section which is reduced compared to the hadron-nucleon cross 
section, the model showed rather good agreement with the available HERMES 
data [T2|13|14] for pions and kaons. In the LUND model, the hadron is formed 
at the formation length Ih 

lh = lp + Zh- ■ (1) 



Here, Ip denotes the prehadron formation or production length. For relativistic 
quarks confined in one dimension the only scale setting parameter is the string 
tension n. Therefore, in the conventional estimate of the Lund model for the 
formation length enters the photon energy divided by the string tension 

. z/ GeV 

L = - , « = !-—. 2 
K fm 



We have fitted the prehadron cross section to the pion-data for the multiplicity 
ratios as a function of Zh and v and found an optimal fit with a prehadron cross 
section equal to (2/3) of the hadron cross section in the extended modelling 
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of Ref. [5] 



(^prehadron ~ g ^hadron- (3) 



This value of the prehadronic cross section is in agreement with Ref. [1] . 

Such a reduction may be partially because the prehadron does not yet have 
the full size of the hadron and therefore interacts with a smaller cross section 
due to color transparency. Besides, a considerable reduction should be also 
expected, since the exponential attenuation of (pre)hadrons used in the fit 
(see below Eq. ([6])) misses the Gribov inelastic corrections [15], which make 
nuclei much more transparent p^l7fl8] . 

Further on, we have considered in [S] the production process for particles 
which cannot be formed from a valence quark in the proton which is knocked 
out by the photon. For example, negative kaons as well as antiprotons can- 
not be formed by a struck valence quark picking up an antiquark from the 
string break-up. They can only be formed from struck sea quarks, which are 
subdominant at HERMES, or from qq pairs formed inside the colour string. 
This different mechanism implies a flavour-dependent formation time from the 
LUND string fragmentation model. 



2 Hadronic p^-broadening as a function of Zh and u 



In this paper we want to test the three stage model further by concentrating 
on quark propagation through the nucleus where most of the transverse mo- 
mentum is acquired. The theoretical calculation is very similar to reference [3]. 
The length Ip after which the prehadron is formed pi5] depends on the energy 
u transferred to the quark, the string tension k and the energy fraction Zh 
of the produced hadron. From energy conservation already follows that if the 
hadron has a very large Zh the quark cannot have radiated very much energy. 
Therefore the formation length of the colour neutral state must have been very 
short [T]|19|20j . The prehadron formation length is computed analytically in 
the framework of the LUND string fragmentation model. If one assumes that 
the prehadrons can be formed directly from the struck quark by picking up 
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an antiquark from the string break up, then the prehadron formation length 
reads 



K 
X 



Zh{l - Zh[ 



Zh 



D. 



2 + Z)a 



iFA2 



1 % 
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Zh 



Zh 



(4) 



Here the parameter Dq is equal to Dq = 0.3 and 2-^1 is the Gauss hyper- 
geometric function. The corrections to the simple Zh{l — 2/i)-behavior of the 
prehadronic formation length Ip given by the Gauss hypergeometric function 
can be recasted into effective powers of Zh and 1 — Zh normalized by an ap- 
propriate prefactor. One can obtain an excellent fit to the "exact" prehadron 
formation length by 



hp — 



1.19 -z 



0.61 



Zh 
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(5) 



We use the length of the quark trajectory Eq. (jlj) to calculate with the dipole 
model the acquired Ap^ of the quark under the constraint that the subsequent 
prehadron is not absorbed on its a way through the nucleus, i.e. that it can be 
finally detected as a hadron. This is necessary, since the information about the 
acquired transversal momentum of the struck quark is encoded in the detected 
hadrons only. 



00 

{Apl)q = {aql)j— J d%dzpA{b,z) J dz'pA{b,z') 

— 00 

• exp —a^ / dz'^ pA (b, 

00 

(S"*) = J d^hdz Pa (b, z^ exp 




(6) 



In this equation, the quantity (crq^) is the mean transverse momentum squared 
qj_ acquired by the quark in one collision multiplied with the corresponding 
quark nucleon cross section. This quantity is related to the dipole nucleon 
cross-section [2T|22] as shown below. Furthermore, we assume a sharp dis- 
tribution of prehadron formation points, namely the prehadron is produced 
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after travelling a distance Ip through the nucleus. Hence, the final induced mo- 
mentum broadening calculated in Eq. can be read as the mean transverse 
momentum q'j_ acquired by the ejected quark in a single collision multiplied 
with the average number of collisions in the nucleus. The resulting transverse 
momentum is averaged over all virtual photon interaction points and weighted 
by the prehadron survival probability. To be precise, the first integral over the 
nuclear density pa averages over all primary interaction points in which a 
quark is ejected from a nucleon. The second integral multiplied with the cross 
section a yields the number of collisions suffered by the ejected quark. The 
exponential factor at the end represents the prehadron survival probability 
S** [5]. It is dictated by the longitudinal thickness of the nucleus at a given 
impact parameter b and by the mean free path of the prehadron in the nu- 
cleus A^^ = cr^ Pa, where the prehadron nucleon cross-section a* = 2/3cr7rAr 
(see Ref. [5]). The mean free path A* is of the same magnitude as Ip and Ra- 
Hence, more weight is given to production points close to the back-surface 
of the nucleus which have large prehadron survival probabilities. In order to 
normalize this expression to the actual number of detected hadrons, we divide 
by the ^-integrated prehadron survival factor (S^,). 

The mean transverse momentum squared times the cross-section, i.e. (crq^), 
can be derived from the dipole nucleon cross section as follows. In the eikonal 
approximation, the ejected high momentum parton moves on a classical tra- 
jectory with impact parameter b and picks up a non-abelian phase factor V{b) 
in the background gauge field generated by the nucleon 



V(b)=V exp 



1^ 



(7) 



Here V{b) is the Wilson line of the parton with impact parameter b relative 
to the proton. We use the notation = A^^t"", where t"'s are the generators 
of the group SU(A'c) in the fundamental representation. The differential cross 
section to produce a parton with transverse momentum q± is given by pro- 
jecting the eikonal phase onto q± and by taking the modulus of the amplitude 
integrated over all possible impact parameters 



da 



d^± (2 it) 



d%d%'e"^^^'-'''^ — (Tr \v\b') V{b) 



N, 



(8) 
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Hence, a fake dipole of size f± = b — b'is constructed from the ejected parton 
in the l^-amphtude and in the l^^-amphtude. Their trajectories are displaced 
from each other by the distance r±. The expectation values of the Wilson lines 
have to be evaluated with respect to the target ground state. In the dipole 
model, the total cross-section for the interaction of a dipole of size •rj_ with a 
target nucleon is given by 



d'b 1- 



Tr 



V\b + f^)V{b) 



(9) 



We define the quantity {crql} as the integral over transverse momentum d'^q± of 
the differential cross section given in Eq. ([8]) multiplied by g^. Differentiating 
the phase factor appearing in Eq. (jH]) twice with respect to the transversal 
separation and performing the integral over (fq± one sees that {crq^} is related 
to the dipole nucleon cross section. 



(2ir) 
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Tr \v\b 



rx=0 



r±) V{b) 

(10) 



This expression confirms the result derived in [21] . Because of the g^-integration 
only the second derivative of the = 0-part in the dipole cross section is rel- 
evant for p_L-broadening. This derivative is a constant due to the color trans- 
parency behavior of the dipole cross section <JdN{r±)rj^-,o oc jl6|. We use 
a form of the x-dependent phenomenological dipole nucleon cross section of 
Ref. [23], which has been adjusted to include soft interactions in Ref. 
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where ro(s) = 0.88 (s/so)""-^^ fm, sq = 1000 GeV^ and 
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Hence, one determines the parameter {crq'j) ^ 4.6 for ^/s ~ 5 GeV, which is 
the typical i/s for the quark nucleon scattering at HERMES energies. 
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Notice that the quark- nucleon differential cross section a is infrared divergent. 
For this reason the mean momentum transfer squared = 0. Nevertheless, 
the broadening, {crq'^), is nonzero and finite. It results from infinitely many soft 
rescatterings with vanishingly small momentum transfers. One can introduce 
an ad hoc infra-red cutoff and regularize the problem. However, this cutoff does 
not affect the final result Eq. flTUl) [21j . The new scale controlling broadening 
is called saturation momentum Qg. For a quark propagating a path- length L 
in nuclear medium of density po the saturation momentum reads 

iQfr = lpoLaoiQfr, (13) 

where (Tq and [Q^)"^ = A/r^ are defined in terms of the saturated form [23|2l] of 
the dipole cross section, Eq. (fTTl) . This increase in saturation scale is identical 
to the broadening of the mean transverse momentum squared of the quark. 
One should mention that Qf is the saturation scale for quarks. For gluons, 
the saturation scale squared is 9/4 times larger. A review of other approaches 
to nuclear broadening can be found in Ref. [25]. In this reference, the relevant 
quantity is qp = {o'q'j_)po which is given as qp = 0.035 GeV^/fm compared 
with our determination qp = 0.032 GeV^/fm at this dipole energy. 

The acquired transverse momentum of the quark given in Eq. ([6]) can be 
computed analytically if one uses a hard sphere approximation for the target 
nucleus which has a homogeneous nuclear density po = (4vr/3rQ) ^ with tq = 
1.2 fm 

PA{b,z)=poQ{RA-b)e{R{b)-\z\) , R{b) = ^Rl - . (14) 

Cylindrical coordinates {b, z) are favorable due to rotational invariance in the 
impact parameter plane. Here, the impact parameter of the initial virtual pho- 
ton is b and Ra = tq ^^^^ denotes the nuclear radius. The hard sphere approx- 
imation is reasonable for large nuclei in which the thickness of the boundary of 
the nucleus is small in comparison to its extension. This approximation gives 
us some insights into the underlying physics of the p± broadening 
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{^pl), = {aql)pMp 1 



{S,)'\SRa QA\Ra 




Q{2Ra-Ip) 



(15) 



The p5_"t>roadening is given by the acquired mean momentum squared per 
unit length (ag^) po times the in-medium propagation length of the quark 
multiplied with the normalized survival probability of the prehadron. This 
product is represented by the expression in square brackets. 

The in-medium propagation length of the quark has two contributions which 
differ depending on the relation between the prehadron formation length Ip and 
the nuclear diameter 2 Ra- If the prehadron formation length is larger than 
the nuclear diameter, then the prehadron is formed outside of the nucleus 
and the in-medium propagation length of the quark is given by the average 
nuclear thickness 3/4 i?^ seen by the quark. The prehadron survival proba- 
bility is equal to one in this case. For prehadron formation lengths which are 
smaller than the nuclear diameter, the in-medium propagation length of the 
quark equals the prehadron formation length Ip plus some higher order correc- 
tions in Ip/ Ra which are due to the finite size of the nucleus. These corrections 
account for the possibility that the prehadron is formed outside of the nucleus 
such that not the entire prehadron formation length Ip contributes. The finite 
size corrections are normalized by the prehadron survival probability, whereas 
the expression which one would expect for an infinitely extended cold nuclear 
medium (the term with Ip) remains unchanged. One should remark, that the 
corrections due to the finite survival probability of the prehadron are numer- 
ically very small in general. If one assumes big homogeneous nuclei for which 
Ra ^ ^p; the finite size effects and the survival probability become negligible 
and the acquired Ap\ is given by 



A hadron with momentum fraction Zh has a {Ap\)h reduced by compared 
to the quark {Ap']_)g. This is a purely kinematical factor and accounts for the 
fact that the average {Apj_)g is shared among the produced hadrons according 
to their energy fractions (c.f. also Ref. [26]). 




{^p±)q = {(^qD^p Po- 
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Fig. 1. p_L-broadening as a function of Zh for pions in Ne, Kr and Xe. 

For a more realistic computation of the transversal momentum broadening of 
the quark, we use a Woods-Saxon distribution for the nuclear density in the 
following. In Fig. [H we compare our results for broadening with the HER- 
MES preliminary data [TT] for vr"*" and vr" for the three nuclei ^°Ne, ^^Kr and 
^^^Xe as a function of Zh- Since the prehadron formation length entering the 
calculation is a function of Zh and we take for the value of u the experimental 
average in the given Zh-hin. As one can see in the plots for the dependence on 
Zh there is qualitative agreement between the calculation and the preliminary 
experimental data. The general shape of the dependence has deficiences: 
For intermediate Zh the agreement is good in all three cases, but in the small 
Zh bin the theoretical Ap\_ for Kr and Xe is too small and in the large Zh bins 
it is too large. Furthermore, our model does not differentiate between vr"*" and 
7T~. In order to do so, one would need to employ a more sophisticated model 
which allows for flavor dependent prehadron formation lengths (see Ref. [5]). 
Within the error bars, the preliminary data do not discriminate between vr"*" 
and 7i~ and the expected effect seems to be small. 

In Fig. [2] we display the dependence of Ap5_ on the photon energy u. Similar 
to the Zh plot, we use in the computation of the prehadron formation length 



10 



the experimental mean value of Zh in the given i/-bin. Because of the increase 
of formation time with u one would expect that the p_|_-broadening increases 
with photon energy. In the prehminary data and in the calculation, however, 
the broadening Ap\ decreases with u. We think that this is due to the ex- 
perimental constraints on the kinematics. With increasing u the experimental 
(zh) (typically z^ ^ 0.35 — 0.45 here) decreases which lowers the resulting 
Ap^. We remark that Ap\ at constant Zh increases with u in the preliminary 
CLAS data at J-Lab [27]. The much lower energy of DIS in this experiment, 
however, may also allow hadron elastic scattering as an important source of 
nuclear broadening. 
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Fig. 2. px-broadening of pions in Ne, Kr and Xe as a function of the photon energy 
I'. 



3 Hadronic p^-broadening as function of 



The variation of p±- broadening with the photon virtuality is the third piece of 
information we have from the preliminary data. Recently an evolution equa- 
tion has been constructed which includes not only the splitting terms in the 
evolution of the parton distributions but also a scattering term for the in- 
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teractions of the parton with the hot medium, i.e. the quark-gluon plasma 
pS] . We will use the equivalent equation in cold nuclear matter. Soft collisions 
with the nucleons do not change the virtuality of the parton, only increase its 
transverse momentum. The "fragmentation functions" D^{zh,Q'^,'P±) give the 
probability for an initial quark q to convert into hadron h with momentum 
fraction Zh, virtuality and transverse momentum p±_ in the course of the 
cascade. The evolution equation for this multiple differential function has the 
following form in vacuum [29] : 



dD^^{zh,Q',p±) a^m }dy 



aiog(Q2) 2n J y 



Pl{y.as{Q')) 



X 



I - y)Q- - f - ,i) A" (f , «^ Px - 

The above equation takes care of the mass constraint y{l — y)Q'^ = Ql/4: + q\ 
arising in the splitting with momentum fractions y and 1 — The transverse 
momenta appear together with longitudinal momentum fractions to guarantee 
boost invariance. After integration over transverse momentum one obtains 
the standard DGLAP equation. For electron-nucleus collisions we consider 
the shower inside nuclear matter with a homogeneous nuclear density po- We 
assume that the medium nucleons change the transverse momentum of the 
quark by giving gj_ kicks, but they do not change its mass scale or virtuality. 
Strictly speaking, this is only true for small momentum transfers i.e. small 
angle scattering. In nuclear matter radiation is interleaved with scattering. 
Therefore a scattering term S{zhiQ'^ ■,p±) has to be added on the right-hand 
side of Eq. (fT8|) . It has two parts, a gain term for scattering into the given 
Zh-hm under consideration and a loss term. 



± 

S{zh,Q^,p±) = ^po J dw J d'^q± 



2h 



da 
d'^q± 



g^px - wq^) - D^^{zn, Q',p±)) 6 



W - Zh 



2mf,z/ I 



(19) 



The evolution allows to calculate two different (p^) from the respective frag- 
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mentation functions in the nucleus and in the vacuum Eq. (fT8|) 



/ d^p^plD{z,Q\p^) 
(pi) = ^4 • (20) 



However, the {p\) defined in this equation is solely coming from the evolution. 
The medium modification of the DGLAP evolution gives the difference of 
mean transverse momentum generated in the evolution in the nucleus and in 
the vacuum. This piece is an additional contribution to the multiple scattering 
contribution {Ap\)h{Q'^) which we fix at = 2.5 GeV^ to the preliminary 
data. In the calculation of the mean transverse momentum broadening of the 
hadron the same averaged (crg^) of the quark and the factor z"^ converting 
quark to hadron transversal momentum squared appears naturally together 
with transverse momentum integrated fragmentation function of this specific 
hadron divided by this hadron multiplicity. 

{AplUQ') = {AplUQ') + zlupo {aql) - ^) • (21) 



To lowest order, Ap\ is generated by the scattering term which gives a higher 
twist contribution to the evolution from to Q^. For > Q^, the evolution 
enhances the mean p\ and for < the devolution decreases the mean 
Although HERMES preliminary data are for relatively small photon vir- 
tualities of = 1.5 — 4.5 GeV^, this yields a sizeable effect for Ap\. As one 
sees in Fig. |3] the calculated Q^-dependence is in good qualitative agreement 
with the preliminary data. It is very encouraging to see that the formalism of 
modified evolution equations which we proposed for the quark-gluon plasma 
can also be related to deep inelastic scattering on nuclei. 

In the plasma, the corresponding scattering term was causing a suppression of 
the fragmentation function in the medium i.e. jet quenching [28]. However, the 
situation in cold nuclear matter is new: The scattering term is much smaller 
and the medium fragmentation functions D{z,Q'^) are almost unchanged. To 
lowest order, the medium- induced (Ap5_)/i(Q^) is originating only from the 
higher twist scattering term oc 1/Q^. Scaling violations due to higher twist 
effects have been considered in [5U|^ for processes with a large Zh- 
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Fig. 3. p_L -broadening of pions in Ne, Kr and Xe as a function of the photon virtuahty 

The dependence of Ap^ on is also discussed in Ref . [32] • This reference also 
suggests two additional mechanisms as possible sources for the Q^- dependence, 
namely NLO processes like photon-gluon fusion and possible colored pre- 
hadrons which lose energy by gluon bremsstrahlung. 



4 Discussion 



We have calculated p^-broadening in transverse momentum distributions from 
the dipole model and compared it to recent HERMES preliminary data. We 
find qualitative agreement with the Zh-, v- and Q^-dependences of p^-broadening 
The dependence on the photon virtuality has been calculated with a modified 
DGLAP evolution equation. Finally we have estimated the effect of absorption 
of the prehadronic state for p^-broadening. 

In a recent paper [32] the nuclear multiplicity ratio Rm has been related to 
broadening, using a similar picture of hadronization [5|5] . In this paper, 
the prehadron formation time is extracted from the multiplicity ratio to tp = 
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Ip oc 0.8^2°'^(1 — Zh)- This formation time is similar to Eq. (jl]) but up to 30% 
smaller at mid Zh- On the other hand the difference between hadronic and 
partonic Ap\ is not spelled out. 

There is another question which needs to be adressed: How do the HERMES 
preliminary data [11] match with the preliminary CLAS data [27]? The main 
difference between CLAS and HERMES is the beam energy, which is 2—5 GeV 
in CLAS in contrast to 7 — 23 GeV at HERMES. Therefore one expects that 
the prehadron formation time Ip is smaller by a factor > 3 at CLAS and conse- 
quently also the resulting hadronic broadening Ap^ would be much smaller. In 
the CLAS experiment, however, effects for hadronic Ap\ are of similar mag- 
nitude as in the HERMES experiment. A possible explanation can be that the 
prehadron stage contributes to the hadronic broadening. At these low energies 
the pion-nucleon elastic cross section is of the same magnitude as the inelastic 
cross section. Therefore elastic scattering competes with absorption for the 
outgoing prehadron. Since the angular distribution of pion-nucleon scatter- 
ing has still sizeable contributions from u-channel exchange, large transverse 
momentum exchanges are possible. A good check is possible when the whole 
angular distribution of the produced hadron is measured. Another important 
feature of the preliminary CLAS data is the linear rise with u which possibly 
saturates at z/ ~ 4 GeV. This linear rise is consistent with Ap\ oc /p oc z/ as 
proposed in Eq. (fT6l) . As discussed, at CLAS hadronization may set in inside 
the nucleus in contrast to HERMES. This does require also a careful analysis 
of the energy dependence of elastic prehadronic scatterings. Therefore in our 
three stage model of hadronization, the second and third step play a more 
important role. 
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